Classification of zero-energy resonances by dissociation of Feshbach molecules 
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We study the dissociation of Feshbach molecules by a magnetic field sweep across a zero-energy resonance. 
In the limit of an instantaneous magnetic field change, the distribution of atomic kinetic energy can have a peak 
indicating dominance of the molecular closed-channel spin configuration over the entrance channel. The extent 
of this dominance influences physical properties such as stability with respect to collisions, and so the readily 
measurable presence or absence of the corresponding peak provides a practical method of classifying zero- 
energy resonances. Currently achievable ramp speeds, e.g. those demonstrated by Diirr et al. [Phys. Rev. A 
70, 031601 (2005)], are fast enough to provide magnetic field changes that may be interpreted as instantaneous. 
We study the transition from sudden magnetic field changes to asymptotically wide, linear ramps. In the latter 
limit, the predicted form of the atomic kinetic energy distribution is independent of the specific implementation 
of the two-body physics, provided that the near-resonant scattering properties are properly accounted for. 

PACS numbers: 03.75.-b, 34.50.-s, 33.15.Fm 



I. INTRODUCTION 

Cold molecules have been the subject of intense study over 
the last few years IliyJ]. In particular, the manipulation of a 
magnetic field around a zero-energy resonance has been used 
for producing both thermal and condensed molecular gases 
from cold atoms. Typically, the magnetic field is linearly 
ramped across the resonance 0, B lj> Li DU V 
alt hough o ther techniques have been demonstrated 
fl^flafl7lfl8ll . Detection of cold molecules often relies upon 
their dissociation, followed by measurement of the resulting 
atomic gas. Dissociation is commonly performed by pho- 
todissociation j^0,[3or linearly ramping the magnetic field 
across the resonance 13, 

Si SI IE mill Ed. 

Linear ramps for molecule detection are often slow in or- 
der to minimise the speed of the fragments to be imaged. 
Such sweeps have also been used to determine resonance 
widths lllll . Fast sweeps have been used to measure the pop- 
ulation of more than one partial wave fl^.fl^l . 

In this paper we study in detail dissociation by linear ramps 
across zero-energy resonances over the whole range of possi- 
ble ramp speeds. We show that the dissociation process can be 
used to probe the nature of the resonance. In the limit of a sud- 
den change in the magnetic field strength the distribution of 
atomic kinetic energy, referred to as the dissociation spectrum, 
can have a sharp peak at an energy dependent on the final 
magnetic field 1 1 1 ] . Our studies show that such spectra occur 
only for resonances with weak interchannel coupling. Conse- 
quently, resonances may be classified by measuring the peak, 
which can be done while little is known about the resonance 
other than its position. This classification scheme applies to 
Feshbach molecules formed from both Bose and Fermi gases. 
The strength of the interchannel coupling affects the suitabil- 
ity of the resonance for studying phenomena of current exper- 
imental interest. For example, halo molecules such as long- 
range dimers and three-body Thomas-Efimov states are more 
easily studied usin g strongly coup led, entrance-channel dom- 
inated resonances Clllalllllllllllll, whereas Feshbach 
molecules created with closed-channel dominated resonances 



are usually unstable with respect to collisions EIH. We 
study the transition from sudden magnetic field changes, re- 
ferred to as jumps, to asymptotically wide, linear ramps, as the 
ramp speed and initial and final fields are changed. We pro- 
vide an analytic treatment of the time evolution for arbitrary 
ramp speeds and initial and final fields. Taking the limit of an 
asymptotic ramp produces a spectrum which has been studied 
both experimentally EEH and theoretically EE El 13. 
We show that two-channel and single-channel approaches pro- 
vide the same asymptotic spectrum. 

In Sections ||T| and HU we introduce the two-channel ap- 
proach to zero-energy resonances and its application to disso- 
ciation spectra, respectively. Dissociation spectra produced by 
magnetic field jumps and their application to classifying res- 
onances are discussed in Section lTVl Section[V]is concerned 
with the transition from jumps to asymptotic ramps. We dis- 
cuss the different theoretical techniques that may be applied, 
and the validity of approximating a given, finite ramp to be 
asymptotic. We conclude in Section lVTl In the appendices we 
present the determinations of the time evolution in the two- 
channel (Appendix lAl and single-channel (Appendix IBl ap- 
proaches. 



II. TWO-CHANNEL SINGLE-RESONANCE APPROACH 

Atom-atom scattering in a cold gas can be magnetically 
tuned because of the coupling that exists between scatter- 
ing channels with different spin configurations. The diatomic 
spin configuration corresponding to the Zeeman state of each 
atom in the gas is referred to as the entrance channel. This is 
coupled to the closed channel, so called because its dissocia- 
tion threshold is far above the relative kinetic energy of atom 
pairs. The coupled-channels character is reflected in the gen- 
eral form of the two-channel, two-body Hamiltonian of the 
relative motion, 



H 2 b(B) = |bg>H bg <bg| + W|bg)<cl| 
+ |clXbg|W+|cl>i/ cl (fi)<cl 



(1) 
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Here "bg" and "cl" indicate the entrance-channel and closed- 
channel spin configurations, respectively. H bg and H c j(B) are 
the Hamiltonians of the two spin configurations in the hypo- 
thetical absence of the interchannel coupling W. We choose 
the zero of energy to coincide at each magnetic field B with 
the dissociation threshold of the entrance channel. Due to a 
difference in magnetic moment between the two channels, the 
closed-channel dissociation threshold can be tuned with re- 
spect to zero energy. Consequently, H C \(B) contains all of the 
magnetic-field dependence of H2r{B). 

Within the magnetic field range of experimental relevance, 
typically only one closed-channel state, referred to as the 
Feshbach resonance state res , has an energy that can become 
near resonant with pairs of atoms in the entrance channel. 
Consequently, the closed-channel Hamiltonian H C \(B) is well 
described by taking into account only res , which fulfils the 
Schrodinger equation 



ft 



(2) 



Here, E les (B) — yU res (fi - fires) is the energy of the Feshbach 
resonance state. The difference in magnetic moment between 
the two channels is denoted yu les , and B res is the magnetic 
field strength where the Feshbach resonance state crosses the 
entrance-channel dissociation threshold, i.e. E KS (B KS ) = 0. 
The single-resonance approach to H C \(B) consists in the ap- 
proximation 

ffci(B) = \^ es }E res (B)(4> KS \ . (3) 
The two-channel Hamiltonian of Eq. ([Q then becomes 



H 2 B(B) = |bg>ff b g<bg| + Wres,bg>^ res ,cl| 

+ |0 res , Cl)<0 res , bg| W + |0 res , cl>£ res (B)<0 res , cl 



(4) 



A zero-crossing of E KS is accompanied by a singularity in 
the scattering length a referred to as a zero-energy resonance, 
at magnetic field Bq. This is shifted from Z? res by the cou- 
pling W. On the side of Bq where a > 0, the highest excited 
vibrational bound state <p b (B) supported by Hre,(B) is usually 
referred to as the Feshbach molecule, and satisfies 



H 2B (B)<p b (B) = E b {E)4> h {B) . 



(5) 



Here E b (B) is the magnetic-field dependent binding energy, 
which vanishes at Bq. In the single-resonance approach, the 
Feshbach molecule assumes the two-component form [26] 



t>2 



(B) 



1 



N b {B) 



G bg (E b (B))W<Pn 



(6) 



Here cf> h h g (B) and cf> c h \B) are the entrance-channel and closed- 
channel components of <p b {B), respectively, and the normali- 
sation constant N b {B) is given by 



(B)= y[ 



1 +• 



WGUE b {B))W\4>^) . 



(7) 



In Eqs. (|6} and Q the energy dependent Green's func- 
tion Gbg(z), associated with the entrance-channel Hamiltonian 



H bg , is evaluated at the binding energy z = E b (B). The so- 
lution of Eq. using the single-resonance Hamiltonian of 
Eq. (0} indicates that E b (B) fulfils the condition 

E b (B) = £ res (B) + <«MWG bg (£ b (fi))W res >. (8) 

The continuum part of the spectrum of HrsiB) is repre- 
sented by the scattering states <p p (B), which are labelled by 
the relative momentum p of an atom pair and obey 



H 2 B(B)(f> p (B) = E4> p (B), 



(9) 



where E = p 2 /m, and m is the atomic mass. The scattering 
states have the following two-component form in the single- 
resonance approach |26|: 



0p E (B) 



+ A(B, E)G bg (E + iO)W<t>« 



(10) 



Here, E + iO indicates that the positive physical energy E is 
approached from the upper half of the complex plane, and 
</>p +) is the background scattering state of relative momen- 
tum p, associated with the entrance-channel Hamiltonian H bg 
via H bg (f> p +) = E0 p + \ The magnetic-field and energy dependent 



amplitude A(B, E) is given by 



A(B,E) : 



<0resF|^ +) > 



E - EM - [A(E) - iy(E)] ' 
where the real functions A(E) and y(E) satisfy 

ME) - iy(E) = (<l> KS \WG b JE + iO)W\<p KS ) . 



(11) 



(12) 



The function A(E) acquires the meaning of an energy shift, 
while y{E) captures the decay width of the metastable Fesh- 
bach resonance state. In the case of negative energies, as in 
Eq. (|8}, the decay width y(E) vanishes. 



III. DISSOCIATION SPECTRA 

Feshbach molecules can be dissociated by changing the 
magnetic field strength to the side of the zero-energy reso- 
nance where the Feshbach molecular state does not exist. We 
consider linear magnetic-field ramps of the form 



B(t) = B as + B(t-t m ). 



(13) 



Here B is the ramp speed, and f res is the time satisfying 
B(t res ) = B KS . The ramp begins at the time t\ and ends at tf, 
ranging between the magnetic fields B(t\) = B\ and B(t{) = Bf. 
The ramp of Eq. Jl 31 corresponds to a linear time dependence 
of the resonance energy 



Eres(t) — E Tes (t ?res) t 



(14) 



where Z? res = yU res B. The validity of Eq. dl4> presupposes 
that the electronic degrees of freedom of the atoms adiabat- 
ically follow the changing magnetic field strength through- 
out the dissociation process. The effect of the ramp on the 
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initial molecular state is described by the evolution operator 
U2B(tf, h)i associated with H^siBit)) by the Schrodinger equa- 
tion 

d 

iti-U 2B (t, t') = H 2B (B(t))U 2B (t, f) , (15) 
at 

with U2 B (t,t) = 1. 

In typical experiments flolfTlll the atomic gas produced by 
the sweep is then allowed to freely expand with the magnetic- 
field strength held at Bf, after which a measurement of the 
atomic momentum distribution is performed. It is therefore in- 
structive to analyse the diatomic state produced by the ramp in 
terms of the scattering states of Eq. evaluated at Bf, which 
then become stationary states. The transitions into these states 
due to the ramp determine the dissociation spectrum n(E), in 
accordance with the general formula | 2(J 

n(E)dE = p 2 dp J </Q p |r dlss (p)| 2 . (16) 

Here J c/Q p denotes the integration over angles, E = p 2 /m, 
and rdi S s(p) is the transition amplitude, given by 

r<Hss(p) = <<|C/2Bfe*i)l4>- (17) 

Here, <f> f p - <p p (Bf) is the scattering state of momentum p at the 
final magnetic field, and <p l h = <pb(Bi) is the Feshbach molec- 
ular state at the initial magnetic field. In the following, we 
consider only resonances for which a spherically symmetric 
Feshbach resonance state is coupled to the entrance channel 
by spin exchange 111 . Consequently, the transition amplitude 
depends only on the modulus of the momentum. 

IV. DISSOCIATION DUE TO A MAGNETIC FIELD JUMP 

In this section we discuss the limit of an infinitely fast ramp 
across a zero-energy resonance, corresponding to an instanta- 
neous change in the magnetic-field strength. This scenario is 
directly related to the recent experiment of Diirr et al. II ill , 
where a sharply peaked dissociation spectrum was observed 
following a fast ramp across the 685 G resonance of 87 Rb. For 
a jump, the transition amplitude of Eq. fl!7i is completely de- 
termined by the overlap of the initial Feshbach molecular state 
and the final scattering state, 

7diss(/>)^<<*pl0b>- ( 18 > 

Evaluating this overlap using the approach of Section [n] re- 
veals that the transition probability density \T& ss (p)\ 2 can have 
a two-peaked structure. An example of such a structure, for 
the 1007.4 G resonance of 87 Rb, is shown in Fig.^ This par- 
ticular form of |7"diss(/?)l 2 can be understood by considering the 
entrance-channel and closed-channel components of 0! and 
<p f p , given by Eqs. and fllOi respectively. The maximum 
at zero momentum is created by the overlap of the entrance- 
channel components of the initial molecular and final scat- 
tering states. The appearance of the peak at nonzero ener- 
gies is controlled by the amplitude A(Bf , E), which determines 



the final scattering states </>L through Eq. (lilt . Provided that 
the influence of A{E) and y{E) in Eq. dl It can be neglected, 
|7diss(p)| 2 tends to peak around E a E res (Bf). This is the case 
for the example of Fig.^ 




P ^ /a Bo J 

FIG. 1: (Color online) Transition probability density |r diss (p)| 2 for 
the dissociation of a 87 Rb2 Feshbach molecule into a pair of atoms 
with relative momentum p, due to a jump across the 1007.4 G res- 
onance. The initial and final fields used here are B, = 1007 G and 
Bf = 1009.4 G. The momentum p mf} = 0.0140 ft/flBohi-, corresponding 
to the resonance energy E mi (B { ) = p^ es /m at the end of the jump, 
gives the approximate position of the peak. 

In general, though, the influence of A(E) and y{E) can be 
significant, leading to the broadening and possibly even disap- 
pearance of the nonzero-energy peak in the dissociation spec- 
trum. For small E = (hk) 2 /m, satisfying \ka\ <K 1, the two- 
body interactions can be completely characterised in terms 
of the scattering length. The regime where the physics de- 
pends only on the scattering length is referred to as the Wigner 
threshold law regime for positive energies, and as the univer- 
sal regime for negative energies. The energy shift A(E) as- 
sumes its maximum value for E — > 0, 

A(E) ~ ^esCBres - B Q ) . (19) 

k—>0 

In this limit A(E) is related to the quantity Bo - B, es , which 
is termed the resonance shift. The decay width y{E) in the 
Wigner threshold law regime assumes the classic form of ll29ll 

y(E) ~ ka bg/ j KS AB, (20) 

k— >0 

where AB is the resonance width, and flb g is the background 
scattering length. The resonance width and background scat- 
tering length are related to the scattering length a by 

a(B) = ahs [ 1 ~B~^B' J- (21) 

Equation (I20> indicates that y{E) can be significant for broad 
resonances with large |flb g |, which will therefore have broad 
peaks in their dissociation spectra. The modulus of the res- 
onance shift Bo - B KS of Eq. dl9t also tends to be large for 
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these resonances. A derivation beyond the scope of this paper 
shows that the resonance shift is given by 1 261 13011 



1 - flbg/a 
1 + fl - «bg/a) 



«b g 

Bo — B KS — AB^— 

a 



Here, a is the mean scattering length given by 1 3 1 
1X3/4) 



a = 



-u 



V2F(5/4) 



(22) 



(23) 



and Z v[ jw = jCwiCg/ft 2 ) 1 ' is the van der Waals length, char- 
acterising the long-range part of the entrance-channel interac- 
tion potential according to V\,Jr) ~ -C(,/r 6 . 

The defining property of entrance-channel dominated reso- 
nances is that the universal formula for the bound state energy, 
£\, ss -H 2 /ma 2 , is accurately extended beyond the universal 
regime by E^, a -h 2 /[m(a - a) 2 ] fil l3lll . The accuracy of this 
extension depends on the degree to which the influence of the 
van der Waals tail of the entrance-channel potential dominates 
over that of the Feshbach resonance state |32]. Conversely, if 
the influence of the Feshbach resonance state is dominant, the 
resonance is referred to as closed-channel dominated. The dis- 
tinction between entrance-channel and closed-channel domi- 
nated resonances is quantified by I20H33I1 



V = 



a tr/(ma ) 



^bg 7-^res 



AB ' 



(24) 



which we term the closed-channel dominance. For entrance- 
channel dominated resonances the relation 77 <K 1 is fulfilled, 
implying that |flb g | and/or |// res AZ?| are large. The latter quantity 
being large indicates strong interchannel coupling 1 26]. This 
leads to the closed-channel admixture of the Feshbach molec- 
ular state being suppressed over a range of magnetic fields that 
is wider than the universal regime. Consequently, the wave- 
function has a long-range halo structure 1 1] within this range 
of fields. 

Although the closed-channel dominance of Eq. ( 1241 was de- 
rived in the context of the properties of the Feshbach molecu- 
lar state, it can be probed experimentally by measuring dis- 
sociation spectra. This can be seen by analysing the posi- 
tion of the finite-energy peak, E pe: ±, relative to its spectral 
width, A£peak. By using the Wigner threshold law expressions 
of Eqs. dl9l and J20i . it is possible to show that the finite- 
energy peak is approximately Lorentzian. Consequently, 
A£peak will be taken as its half-width at half-maximum. A 
distinct peak appears in the jump dissociation spectrum when 
£peak/A£peak » 1. Evaluating this quantity in the Wigner 
threshold law regime yields the expression for the peak clar- 
ity, 



J peak 



PresC^f _ Bq)/E. 



vdW 



l/(2/7 2 ) 



A£peak V/WBf - fl )/£vdW " V(V) 



(25) 



where we refer to £Vdw = ft 2 /ma 2 as the van der Waals en- 
ergy. £ v dw, which greatly exceeds the relative kinetic ener- 
gies typical of cold atom pairs, gives the order of magnitude 



of the spacing between the highest excited vibrational bound 
states. It therefore sets the largest energy scale able to be 
experimentally accessed without other resonance phenomena 
arising, and limits the range of energies able to be considered 
within a single-resonance approach. In the range of validity 
Hks(B{ - Bo) < £\,dw, Eq. fl25i gives the inequality 



^peak 



AE, 



i-iW) 

peak ' yjl - 1/(4t/ 2 ) 



1 



(26) 



with equality for yU res (Bf - Bo) = £Vdw- It is evident from 
Eq. d26i that the condition for the existence of a distinct peak 
at a nonzero energy is 77 » 1, which is characteristic of closed- 
channel dominated resonances. An example of a resonance 
falling into this category is provided by the 1007. 4G reso- 
nance of 87 Rb, for which 77 = 5.9 (see Fig. 0. Entrance- 
channel dominated resonances, characterised by 77 <K 1, ex- 
clude the possibility of a clear peak appearing at a nonzero en- 
ergy in the dissociation spectrum. The presence or absence of 
a sharp peak in the dissociation spectrum is therefore an indi- 
cator of whether a zero-energy resonance is entrance-channel 
or closed-channel dominated. 

We have verified the quality of the estimate given by 
Eq. d25i by comparing it to the results of calculations us- 
ing the approach of Section [H] The comparison, shown in 
Fig. confirms that Eq. fl25l > provides a good estimate of 
the ratio £p e ak/A£p e ak- Moreover, the analytic expression of 
Eq. d25l > slightly underestimates £p e ak/AEp e ak, meaning that 
the peak clarity one would expect to see experimentally is al- 
ways slightly greater than the analytic estimate. 




689 



FIG. 2: (Color online) The measure of peak clarity Epeak/AEpeak 
in the dissociation spectrum of a 87 Rb 2 molecule, using the 685 G 
resonance. This resonance has a high closed-channel dominance of 
77 = 463. The initial field used is B; = 684.5 G. The dots indicate val- 
ues extracted from fits to the spectra obtained from full two-channel 
single-resonance calculations, while the solid line shows the approx- 
imate analytic predictions of Eq. <25> . At B ~ 688 G, the single- 
resonance approach breaks down as the energy range covered starts 
to exceed the van der Waals energy £ v dw • 



The preceding analysis treats the ideal case of an instanta- 
neous change of the magnetic field strength. An experiment 
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is described in Ref. full that actually operates in this regime. 
After an 80 G/ms ramp of the magnetic field across the 685 G 
resonance of 87 Rb, the atomic gas was found to expand out- 
wards as a spherical shell, corresponding to a sharply peaked 
dissociation spectrum. In order to assess how closely the 
jump approximation follows the experimental observations of 
Ref. flUl . we have performed an exact two-channel calcula- 
tion of the dissociation spectrum of Eq. dl6> . for the exper- 
imental ramp speed of 80 G/ms. The comparison is shown 
in Fig. |5] and confirms that the 80 G/ms ramp is fast enough 
to produce a spectrum virtually indistinguishable from that of 
a jump. Consequently, the experiment of Ref. lUUl provides 
a clear probe of the closed-channel admixtures of the initial 
bound and final scattering states, and of the strong closed- 
channel dominance (77 = 463) of this particular resonance. 
We note that experiments reported in Ref. II ill using the less 
closed-channel dominated (77 = 5.9) 1007 G resonance of 87 Rb 
did not produce such a peak for the 80 G/ms ramp-speed used. 
We have calculated that to be in the jump regime for this res- 
onance would require a ramp speed of 5000 G/ms or above. 
The minimum ramp speed to be in the regime of a jump is 
different for each resonance, however ramp speeds of order 
1000 G/ms are now possible |3- 
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FIG. 3: (Color online) Dissociation spectrum of 87 Rb2 molecules 
produced in the vicinity of the Bo = 685 G zero-energy resonance. 
The initial and final magnetic field strengths are Bj = Bq - 50 mG 
|H andfi, = Bn + 40mG ITHl. 80 G/ms is the ramp speed reported 
in Ref. 1 11], while 1000 G/ms is approximately the maximum cur- 
rently available 1 34]. The distinct peak is centred around h x 0.073 
MHz, corresponding to 3.5 fiK, as compared to the 3.3 [iK that can 
be deduced from Ref. llUl . 



of Rb and found it to be negligible. For the 48 G resonance 
of 133 Cs referred to in Fig. 0] however, we have found that 
the transfer into the bound state can be as much as 10 %. In 
our calculations we have assumed sweep rates of no more than 
1 000 G/ms, which is currently achievable experimentally 1 34 ] . 
The transfer is unusually large for this particular resonance 
because the entrance-channel potential supports a state with 
an energy of only £_j = -h x 45 kHz This makes the 
avoided crossing between the bound states narrow, as shown 
in Fig.@] and so the wave functions of the Feshbach molecule 
and the final bound state converge to (/>_i quickly as B moves 
away from Bq. This increases the overlap between the initial 
and final bound states and so makes the transfer from a jump 
larger. 
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FIG. 4: (Color online) Energies of the two highest excited vibra- 
tional bound states near the 48 G resonance of 133 Cs. For fields far 
below the resonance, the Feshbach molecular state adiabatically cor- 
relates with the highest excited bound state </>_[ of the entrance chan- 
nel (dotted line), while the next deeper bound state corresponds to 
the Feshbach resonance state (dashed line). Above B the Feshbach 
molecule does not exist and the next deeper state converges to 
Due to this avoided crossing being narrow, the possible transfer into 
the final bound state from a jump can be as large as 10%. 



TRANSITION FROM JUMPS TO ASYMPTOTIC 
RAMPS 



At magnetic fields far from Bq, the highest excited vibra- 
tional bound state of H^aiB) often corresponds to the high- 
est excited vibrational state of the entrance channel potential, 
<f>-\, as illustrated in Fig.@] When the Feshbach molecular 
state is subjected to a magnetic field sweep across the reso- 
nance, some population can then be transferred to the highest 
excited vibrational bound state at the final field strength. Some 
Feshbach molecules will therefore not be dissociated during a 
ramp. We have calculated this effect for the 685 G resonance 



In this section we consider ramp speeds that are intermedi- 
ate between the jumps of Section HV1 and asymptotic ramps. 
Keeping the initial and final magnetic fields shown in Fig. [3] 
for all ramp speeds, we calculate the dissociation spectra in- 
cluding the full time evolution of Eq. J17I . The resulting spec- 
tra have the characteristic double-peaked structure for jumps, 
but coincide with the asymptotic form for sufficiently slow 
ramps. This transition is shown in Fig. [5] For any given ramp 
speed, however, making the ramp wide enough will result in a 
spectrum of the asymptotic form. 
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FIG. 5: (Color online) A sequence of dissociation spectra of 87 Rb2 
molecules, produced near the 685 G zero-energy resonance, for a 
range of ramp speeds. The initial and final magnetic field strengths 
are B, = B - 50 mG and B f = B n + 40 mG, respectively, and cor- 
respond to the conditions of Ref. 1 1 1]. The dash-dotted curve corre- 
sponds to a jump of the magnetic field strength and is identical to the 
curves shown in Fig. [3] The circles illustrate the spectrum obtained 
from an exact two-channel calculation for B = 0.1 G/ms, and coin- 
cides with the evaluation of Eq. 1321 indicated by the superimposed 
solid line. 



For linear ramps of the form of Eq. dl 31 . it is possi- 
ble to derive an exact formula for the time-evolution op- 
erator UzsitJ 1 ) associated with the two-body Hamiltonian 
HzR(B{t)). Here we consider a continuum of positive ener- 
gies as represented by the scattering states of Eq. I ll PL in con- 
trast to the set of discrete levels examined in previous analytic 
treatments Ill l35ll36ll . We separate Eq. @ into its stationary 
and explicitly time dependent parts, 



H 2B (B(t)) = ff stat + H d {f) , 



(27) 



where i/ stat = H2b(B k& ), and we now express the single- 
resonance, closed-channel Hamiltonian of Eq. Q as a func- 
tion of time, i.e. H c \{B(t)) — > H c \(t). Using the separability 
of H c i(t), the derivation of Appendix 1X1 leading to Eq. ( IA5i 
shows that U2B,(t,t') can be calculated from only its closed- 
channel component, 



1 



*ZB % f '> = fo Q(t ~ Ofa*. C W2B(', Ol&es, cl) . 



(28) 



Here, 6(t - t') is the step function, which yields unity for 
t - f > and elsewhere. The linear variation of the reso- 
nance energy makes it possible to represent (t, ?') in terms 
of the following Fourier integrals: 



C e -iE('-r')/n nco 

£ ,ft0= J^-2^JL dE ' 



2 /[<Kz'/)-^(z/)] 



(29) 



Here z — E + iO, z' — E' + /0, and the phase <p(z, t) is given 
in Eq. jA8> . This leads to the expression of Eq. ( IA11H for the 
time-evolution operator. Consequently, U2B(k, t\) is expressed 



solely in terms of the known energy-dependent Green's func- 
tion G s tat(z) associated with the stationary Hamiltonian H stat 
and its closed channel matrix element g s tat(z)- 

In Appendix|X]we use Eq. JA10I > together with a stationary 
phase analysis to evaluate the transition amplitude of Eq. dl7t 
in the limit of an asymptotically wide ramp. This shows the 
transition probability density to be 



I^ssCp)! 2 = 



In 
x exp 



K^l WreJI 2 X 



2tt 

h\E Ki 



j dp'8(p-p')\(ffiw\(P^)\ 2 ] ■ 



(30) 



Under typical experimental conditions the energies of the dis- 
sociated atoms do not exceed several jjK, which is usually 
well within the Wigner threshold law domain. Consequently, 
the interchannel coupling matrix elements in Eq. J30l > can be 
evaluated in the limit of zero relative momentum and related 
to the resonance width AB through |2€ 



AB 



(2nK) 3 



(+)\ 



47r/i 2 flb e 



(31) 



This constitutes a rigorous derivation of the asymptotic disso- 
ciation spectrum, 



n(E) 



d 

-dE™* 



E\AB\ S 




(32) 



which was previously inferred in Refs. dues from a Fermi 
Golden Rule argument. The accuracy of this expression was 
confirmed by the experiments of Mukaiyama et al. 11011 with 
23 Na and Diirr et al. 1 1 1] with 87 Rb. This led to the applica- 
tion of Eq. d32l > to measuring the widths of previously unex- 
plored resonances II ill . The dissociation spectrum of Eq. J32t 
is plotted in Fig. [5] as the asymptotic limit of slow ramps for 
the 685 G resonance of 87 Rb. 

The dissociation spectrum in Eq. d32l > does not depend on 
the magnetic moment fi Tes of the resonance state, but just on 
|flb g AB|. This product characterises the scattering length in 
the close vicinity of a zero-energy resonance according to 
a » -a\, g AB/(B - Bo). This suggests that it should be possible 
to arrive at Eq. d32i using different descriptions of the underly- 
ing two-body physics, the only condition being that the near- 
resonant scattering properties are correctly recovered. Single- 
channel descriptions 1321 1371 13811 are capable of describing 
entrance-channel dominated resonances, such as the 48 G res- 
onance of 133 Cs shown in Fig. 0J over a wide range of mag- 
netic fields. The approach of Appendix 151 varies the poten- 
tial in such a way that the scattering length a is recovered for 
all magnetic fields, as well as the energy of the high- 
est excited vibrational entrance-channel state away from the 
resonance. For closed-channel dominated resonances, how- 
ever, single-channel approaches are valid only in the largely 
inaccessible universal regime. Taking a separable form of the 
potential allows Eq. d32i to be derived analytically, as shown 
in AppendixlBl Consequently, for asymptotic ramps the phys- 
ical description need only apply to the near resonant region, 



7 



whereas for jumps it must apply over the entire range of mag- 
netic fields. 



VI. CONCLUSIONS 

We have studied in detail the dissociation spectra produced 
by linear sweeps of the magnetic field strength across zero- 
energy resonances, for field ranges and ramp speeds varying 
from jumps to asymptotic ramps. A key finding of our stud- 
ies is that jumps in the magnetic field can be used to classify 
zero-energy resonances according to the dominance of their 
entrance and closed channels. Our scheme requires measuring 
whether a sharp peak occurs in the dissociation spectrum near 
the final resonance energy. The presence or absence of this 
peak allows a resonance to be classified as closed-channel or 
entrance-channel dominated, respectively. This may be mea- 
sured before other properties of the resonance are studied. By 
comparing the spectrum of a jump to exact numerical calcula- 
tions for the example of the 685 G resonance of 87 Rb, we have 
shown that the ramp speed required for such a measurement 
is well below the maximum currently possible. 

Our calculations of the dissociation spectrum for varying 
ramp speeds and fixed initial and final fields show that the 
spectrum approaches the asymptotic form as the ramp speed 
is lowered. Increasing the width of the ramp for a fixed ramp 
speed also brings the spectrum closer to the asymptotic result. 
We have analytically determined the time evolution for linear 
ramps using a two-channel approach. Evaluating this in the 
asymptotic limit rigorously produces the asymptotic spectrum 
of Eq. (I32i . Asymptotic ramps, however, only require the 
physical description to correctly describe the near-resonant re- 
gion, and so may also be described using a single-channel ap- 
proach. 
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APPENDIX A: TWO-CHANNEL TIME EVOLUTION 
OPERATOR 

In this appendix we derive an exact expression for the time 
evolution operator in the two channel, single resonance ap- 
proach. We then use this result to evaluate the transition am- 
plitude for the case of an asymptotic ramp, which provides 
a rigorous derivation of the corresponding dissociation spec- 
trum. We calculate the transition amplitude in terms of the 
retarded Green's function G^(t, t'), related to the time evolu- 



tion operator U2B(t, O by 



G%(t,t') = -6(t-t')U 2 B(t,t'). 



(Al) 



We separate the constant and time-dependent parts of the 
Hamiltonian as in Eq. (127 \ . The retarded Green's function 
associated with the stationary Hamiltonian H stat in Eq. i21\ is 
given by 



G ta t (f - t') = -9(t - t') exp[-/// stat (f - t')lh] . 
in 



(A2) 



The separation of H2b(B) into stationary and time-dependent 
parts then gives two equivalent integral equations for the full 
retarded Green's function: 

G< + B ) (r,(') = O-0 + J dTG^(t,T)H cl (T)G ( ^ t (T-t'), 

(A3) 

referred to as the preform, and 

G«(f, ?') = G«(f -t') + J dTG ( ^(t - t)H c1 (t)G<£(t, ?') , 

(A4) 

termed the postform. Differentiating either form of the in- 
tegral equation with respect to t and using Eq. flAlb recov- 
ers Eq. (TBI defining Uzsit, t'). Substituting Eq. ( IA3I into 
Eq. ( IA4I > and using the form of H c \{t) in Eq. Q gives the fol- 
lowing general expression for the two-body time-dependent 
Green's function: 

G«(ff,fi) = G«(? f -fi) + 



+ </TG^ t (f f -r)|0 res ,cl>£ res (T) 



<fc»,cl|G£i(T-n) + 



■/ 

j dr' 8 %(t, T')£ re s(T'X0 res , c!|G«(t' - ti) 



(A5) 



Here, g^(f, f) is given in Eq. j28l >. Throughout these appen- 
dices we use the notation of a lower-case g to refer to the 
closed-channel matrix element of the corresponding Green's 
function. 

The explicit form of G^ ){tf, tO in Eq. ( IA5> shows that ob- 
taining its closed-channel matrix element g^(t,f) is suffi- 
cient for fully determining the time evolution. To this end, 
we project Eq. JA4> onto |0 res , cl) from the left and right and 
again use the definition of H c \(t) to yield 



,(+) 



,(+) 



f) + j drg^ 



(f-T)E„(T)*&W) 



(A6) 



Performing a Fourier transform of both sides of Eq. JA6t with 
respect to t leads to the time dependence of E KS (f), given in 
Eq. dl4> . being converted into a derivative with respect to en- 
ergy, yielding the following ordinary differential equation: 



in—giaiz, ?') = n 



BE 



<piz, 



S2b(z, t') + l/E n 



(A7) 
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Here the energy and time dependent matrix element g2h(z, t') 
is given by a Fourier transform with respect to time, 
82b(z, f) = fdt e iz( '-'' )/t, g ( ^(t, t') . The argument z = E + ;0 
has the dimension of energy and its imaginary part provides 
the appropriate convergence factor. The energy and time de- 
pendent phase ip(z, f) is given by 



<p(z,t') 



E( f - t KS ) 
h 



1 r E 

frE Tes Jo 



dE' 

gstat(z') 



(A8) 



where z' = E' + iO, and g sm (z) = J dte^'-^g^t - t'). 
A solution of Equation (|^) can be inferred from the deriva- 
tive with respect to energy of e'^' '#2b(z, f). Calculating the 
inverse Fourier transform of the relation for g2B(z, f) thus ob- 
tained gives Eq. J29I . 

Substituting g2B(t, f) into Eq. ( IA5> leaves the resonance en- 
ergy and the phase terms as the only explicitly time dependent 
quantities. The time-dependence due to the resonance energy 
may be eliminated by writing 

(r - f res y Vfcf) = -m^ <M - ) JL e w-t,«)/n (A9) 

oE 

and calculating the partial integrals with respect to energy. 
The time integrals in Eq. dA5l > may then be expressed as 
Fourier transforms, reducing the general expression for the 
two-body retarded Green's function to 



. r dE 



,-iE(ft-fd/h 



G sta t(z) X 
" I0res,cl) 



1 



f 

<J —Cx 



dE'\<p KS ,cl) 



gstat(z) 
il<p(zA)-<P(z' A)] 



<0 res ,cl|G smt (z) 



gstat(z)^-E , resgstat(z') 



s ,cl|G stat (z')j 
(A10) 



This Green's function represents the exact evolution of the 
two-body system in the two-channel, single resonance ap- 
proach. The separable form of the closed-channel potential 
and the linear time dependence of the resonance energy allow 
the dynamics to be calculated from the known Green's func- 
tion G s tat(z) associated with the stationary Hamiltonian // stat 
and its closed-channel component. 

We will now evaluate the transition amplitude T^ip) of 
Eq. fl!7l > for the case of an asymptotic ramp, which allows the 
integrals contained in Eq. JAlOt to be calculated analytically. 
In the limit of fj — > -oo, £ , res (f;) is large and negative, and so 
we can write 



G s tat(z) 



1 - |0 res ,cl> 
G st at(z) 



1 



/- — > — CO 

= G 2B (Bi,z). 



gstat(z) 

1 + |0 res ,cl> 



<0 res ,cl|G sta t(z) 

E ies (tj) ( 
1 - E KS (ti)g stzt (z) 



a, cl|G sta t(z) 



(All) 



The energy-dependent Green's function G2b(B{,z) is associ- 
ated with H2B(Bi). We note that since the singularities of 



G2b(B,z) determine the two-body energy spectrum, the mag- 
netic field dependent energy of the molecular bound state is 
constrained by the condition 



-^res (B)g s UE b (B)) = 1 



(A12) 



which is equivalent to Eq. The contribution given by 
Eq. iAl 1> to the transition amplitude from the initial bound 
state (p l b to the final continuum state ^j, is then 



zM) 



e 



-«W*Xft-*)/»<^|^> 



(A13) 



In the asymptotic limit of t\ — > -oo and tf — > +oo, the in- 
terchannel coupling becomes negligible at the initial and final 
fields. The initial Feshbach molecular state 0? and the final 
continuum state 0j, are then orthogonal, and the contribution 
of Eq. iAl 1> to T& m (p) vanishes. 

The remaining term in Eq. jAlOt contains the rapidly vary- 
ing phase terms <p(z', t\) and ip(z, tf), which we evaluate using 
the stationary phase method 1 1 ] . First we evaluate the matrix 
elements that appear in the integrals over E and E'. The point 
of stationary phase for the term <p(z', h) occurs at z' = E\,(Bi). 
Using the definition of the Feshbach molecular state in Eq. (|6jl, 
and Eqs. and iAXl we find that 



<0 ies , cl|G stat (z')|0b) 

gstat(z') 



N h (Bi) 



= ^-hE Ksl p"(Elti). (A14) 



Here we have also made use of the following representation 
of G stat (z) in terms of the entrance-channel Green's function 
Gbg(z) and the coupling W: 

1 

G st at(z) = |bg)G bg (z)<bg| + |0res, cl> — <0 res , cl| 

z 

+ G bg (z)W|0 res ,bg)<0 res ,cl| + |0 res ,cl>^ res ,bg|WG bg (z) 

z 

+ |0 res , cl)i(0 res |WG bg (z)W|0 res )<0 res , cl|G sta t(z) • (A15) 

z 

Using the definition of the scattering states in Eq. dlOt we also 
find, in the asymptotic limit of tf — > +oo where the interchan- 
nel coupling can be neglected, 



«|G stat (z)|^ es ,cl) 
gstat(z) 



« +) |W|0res) 

t f -»+oc i — p 2 /m 



(A16) 



The results of Eqs. ( IA14I and ( IA16I then give the transition 
amplitude to be 



*(+) 



7diss(» = - 



|W|0 res >A/b(fli) 

: 1 

2nHE KS 



-/ 



dE 



-itp(z,h) 



z - p 2 /m 

dE'e-^'^''^-^ . (All) 

<J — OO 
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Reversing the order of integration, the energy integrals can be 
evaluated in the limit of f; — > -oo, giving 



TdissO?) = - 



2n 



•<< ) |W|0 res >x 



x exp (i[ip{E h {Bi), fj) - if(p 2 /m + iO, f f ) + tt/4]) . 

(A18) 

The imaginary part of the phase difference will give a real 
damping factor when we calculate |7di ss (p)| 2 , and is there- 
fore the quantity of interest. Using the spectral decompo- 
sition of Gbg(z) to evaluate the imaginary part of the phase 
factor in Eq. JA1 81 then gives the result of Eq. (I30i . which 
leads to Eq. fl32i for energies E in the Wigner threshold law 
regime. We note that the accuracy of Eq. i32\ depends on the 
quality of the asymptotic evaluation of the phase integral in 
Eq. jA18> . For a given resonance and initial and final mag- 
netic field strengths, the associated conditions are easier to 
fulfil for slower magnetic field ramps. 



APPENDIX B: ASYMPTOTIC DISSOCIATION SPECTRUM 
IN THE SINGLE-CHANNEL APPROACH 

The effective single-channel approach aims to recover the 
resonance-enhanced scattering properties by choosing the in- 
teraction potential as a perturbed background scattering po- 
tential. If such an interaction potential is taken to have a sep- 
arable form, the two-body Hamiltonian can be chosen to be 

Hi! 

H 2B =Ho + VeffCB) ■ (Bl) 
Here the associated effective interaction potential is 

Veff(B) = U-b g >£(B)trbgl , (B2) 

and Ho = -h 2 V 2 /m is the kinetic energy operator of the rel- 
ative motion of an atom pair. V e ff(B) is characterised by the 
form factor Xbg and the amplitude t;(B). The magnetic-field 
dependence of V e ff(B) is contained in £ {B) through its depen- 
dence on the scattering length a(B) 1 37 ] : 



- Ctbg|Go(0)^ b g> + 



4jrh 2 a(B) 



(B3) 



Here Gq(z) = (z - Hq)~ 1 is the free Green's function associ- 
ated with the Hamiltonian Ho and evaluated at zero energy in 
Eq. JB3> . while Ott>g|0) denotes the overlap of the form factor 
Xbg with a zero-momentum plane wave. 

Given the resonance-enhanced behaviour of the scattering 
length of Eq. fl2 li . Eq. jB3> in general implies a nonlinear de- 
pendence of %(B) on the magnetic field strength B. However, 
if B varies linearly with time, it can be shown that in the close 
vicinity of the singularity of the scattering length at B = Bo, 
the amplitude of the separable potential follows the linear 
dependence given by 



#t)=&+f(f-A>). 



(B4) 



Here 1/fo = (Xb g \Go(0)\xb g ), ?o denotes the time when the 
zero-energy resonance is crossed, and 



mB(27:hy\(x bg \0)\ 
47rh 2 a hg AB\(x hg \Go(0)\Xb g )\ 2 



(B5) 



In the following, we will assume that the linear form of 
Eq. (IB4> applies throughout the whole dissociation ramp. It 
can be shown that the final result is independent of the precise 
time dependence of £(f) outside the near-resonant region. 

Given the linear time dependence of Eq. iBAl and the sep- 
arable form of the interaction potential in Eq. dB2b . the effec- 
tive single-channel Hamiltonian is similar to the two-channel 
Hamiltonian of Eqs. ( I27i and Q. Consequently, the analytic 
determination of the evolution operator for the single-channel 
Hamiltonian of Eq. (IB li proceeds along the lines of Appendix 
|X]by noting the following substitutions: Hq plays the role of 
Hgtut, while £res(0 is replaced by ^(f) and |</> res , cl) by |^bg)- A 
similar substitution is applied to the respective Green's func- 
tions and their matrix elements. 

We will now evaluate the general expression of Eq. dl7> for 
the dissociation transition amplitude in the asymptotic limit 
of f; — > -oo and tf — > +oo for the single-channel Hamilto- 
nian. Applying the above substitutions, the contribution of 
Eq. iAl 1> to 7di ss (p) vanishes for the same reasons as in the 
two-channel case. The first phase integral of the remaining 
term in Eq. is analogous to the two-channel case, and 

is given by Eq. JA14t with the appropriate substitution. The 
second phase integral (corresponding to Eq. ( IA16i in the two- 
channel case) takes the form 



f°° dE 
J Eh{Bi) 2nh' 



<^ +) |G (z)bTbg> 



-Mz,tf) xy P 



go(z) 



~ --e 

f r ->oo Jl 



-i<p(p~ /m+/0,/f) 



<Pkbg> 



go(p 2 /m + iO) 
This leads to the following result for Tdi ss (p): 



(B6) 



TdissO?) 



<plW 



X 



m go(p 2 /m + 20) 
x exp (/[^(£ b (Si), fi) - tp(p 2 /m + iO, t t ) + tt/4]) 



(B7) 

To calculate the imaginary part of the phase difference we 
use the definition of tp{E, t) of Eq. JA8> . yielding 



In#(£ b (Bi), k) - V(p 2 /m + *0, f f )] 



—rim I 

H Jo 



dE 



j(E + iO) 



(B8) 



Multiplying the numerator and denominator of the integrand 
by -£(tf), we then take the asymptotic limit of ff — > +oo, in 
which case £(ff) is large and we can write 



Im[^(E b (Bi), fi) - ip(p 2 lm + iO, f f )] 



1 



*->«> ^go(O) 



-Im 



Jr*p 2 /m 
dE- 
o 1 



l-g(r f )g (0) 
l~£(tf)8o(E + iQ) 



(B9) 



10 



We then make use of the general form of the T-matrix of any 
separable potential lErill . 



<o|r eff (£ + /o)|o> 



g(ff)l(y bg |0)l 2 

l-{(tt)go(E + iO) 



(BIO) 



Together with the explicit expression for £ given in Eq. ( IB5I . 
and the low-energy expansion of the T-matrix, this gives the 
following result in the Wigner threshold law regime: 



Im ^(£ b (Si), k) - <p(p 2 /m + M, f f )] = - 



2 a hg AB 

3 H 2 mB ' 



(Bll) 



By evaluating the remaining prefactors in Eq. dB7l in the zero- 
momentum limit, we arrive at the final formula for the disso- 
ciation transition probability density, 



l^dissC?)! 2 



which recovers Eq. ( 1321 



|flb g Afi| 
nh 2 m\B\ 



exp 



4 flbgAB 
~3 tihnB 1 



(B12) 
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